Electronic structure of the CasCc^Og compound from ab initio local interactions 
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We used fully correlated ab initio calculations to determine the effective parameters of Hubbard 
and t - J models for the thermoelectric misfit compound Ca3Co4C>9. As for the Na x Co02 family the 
Fermi level orbitals are the dig orbitals of the cobalt atoms ; the e' g being always lower in energy by 
more than 240 meV. The electron correlation is found very large U/t ~ 26 as well as the parameters 
fluctuations as a function of the structural modulation. The main consequences are a partial a\ g 
electrons localization and a fluctuation of the in-plane magnetic exchange from AFM to FM. The 
behavior of the Seebeck coefficient as a function of temperature is discussed in view of the ab initio 
results, as well as the 496 K phase transition. 



I. INTRODUCTION 

First studied for their ionic conduction properties^, 
the layered cobalt oxides have been, in the last years, 
the object of a regained attention, due to the discov- 
ery of peculiar transport and magnetic properties. One 
can cite, for instance, the superconducting state of the 
Nao.35Co02 — 1.3H 2 0j2 the large thermoelectric power 
(TEP) in the Nao.7Co02j^ or the rich phases diagram of 
the Na x Co02 family^ In terms of applications the large 
thermoelectric power (TEP) found in some materials of 
this family is certainly the most interesting property. 
Indeed, in addition to Seebeck coefficients larger than 
100/iV/K at room temperature, these systems present 
the atypical association of a low resistivity and a low 
thermal conductivity. It results a large figure of merit 
ZT = S 2 T/pK at room temperature^ (where Z is the 
factor of merit, S the Seebeck coefficient, T the temper- 
ature, p the electrical resistivity, k the thermal conduc- 
tivity) . 

The motivation was at the origin the search for new 
compounds, based on the same C0O2 layers (known 
to be responsible for the electronic and magnetic 
properties), but with better chemical stability than 
the Na x CoC>2 ones. Indeed, in air, these systems can 
both easily accept water or other small molecules 
in between the C0O2 layers and loose part of their 
sodium content, both processes affecting the desired 
properties. The main idea was to replace the alkaline 
layer by a more complex oxide one. Many systems were 
thus synthesized such as Tl a [(Sr, Ca)i_^]i +x Co02r 
Bi Q [Ao. 7 5Bi .25 0]3 + 3 x/2 Co02 (A Ca,Sr)£, 
Sr 2 2 Co02^ etc. 

Among these layered cobalt oxides, the so-called 
CasC^O^ ([Ca2Co03]o.62[Co02]) occupies a special 
place. Indeed, this chemically "simple" system has been 
used as the reference system for the study of the ther- 
moelectric properties in Co02-based compounds. In ad- 
dition to its interesting thermoelectric properties the 
Ca3Co40g present puzzling magnetic transitions^ with 
the inset of an incommensurate spin density wave (IC- 
SDW) at low temperature, an often called "ferrimag- 
netic" transition around 19 K and a transition inter- 



preted as a cobalt spin state transition around 380 K. 
The long-range IC-SDW sets in place at 27 K however 
from 100 K a short range IC-SDW order is observed by 
muon relaxation^. 




Figure 1. (Color online) Structural representation of the 
Ca3Co40g compound. Green : C0O2 layers of edge-sharing 
octahedra, blue : cobalt octahedra in the rock-salt layers and 
light gray : calcium ions. 



The Ca3Co40g system is built from the alternation of 
Cdl2-type C0O2 layers and CoCa203 rock-salt-type lay- 
ers^ (see fig. |TJ). These layers are stacked along the c-axis 
and present incommensurate lattice parameters along 
the b direction. The C0O2 layers are built from edge- 
sharing cobalt octahedra presenting a compression along 
the c crystallographic direction (the octahedra (111) di- 
rection). The rock-salt layers are three- fold layers built 
from two Ca-0 planes sandwiching a Co-0 plane, or 
cquivalently built from a layer of cobalt octahedra, as in 
a rock-salt structure. The incommensurability between 
the b lattice parameters of the two subsystems induces, 
in each layer, a structural modulation presenting the pe- 
riodicity of the other subsystem. One of the important 
consequences is a distortion of the coordination sphere 
of each cobalt atom in the C0O2 layers (different Co-Co 
distances, different Co-0 distances and Co-O-Co angles). 
Indeed, the C0O2 subsystem has been shown to be re- 
sponsible for the low energy properties of the compound 
(transport, TEP, magnetism, etc.), since it supports the 
Fermi level electrons^^ while the CoCa203 subsystem 
was proved to be gapped. In addition, in the related 
Na x Co02 family the Fermi level properties were proved 
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to be strongly dependant both on the doping^ and on the 
local structural parameters^ 3 -. 

The purpose of the present work will thus be to study 
the effect of the incommensurate modulations on the low 
energy degrees of freedom of the CaaC^Og system. At 
this point one should notice that, unfortunately, density 
functional (DFT) electronic structure calculations fail to 
properly describe the electronic structure of the C0O2- 
based systems. Indeed, LDA as well as LDA+U calcula- 
tions^ exhibit a conducting CoCa 2 3 subsystem dom- 
inating the density of states at the Fermi level. DFT 
is known to illy treat strong correlation effects, how- 
ever, in many strongly correlated systems, it is never- 
theless able to yield reasonable Fermi surfaces and the 
correct magnetic ground state. In the CoC^-based sys- 
tems, this is not the case. In the Na x CoC>2 compounds 
DFT calculations yield a ferromagnetic ground stated, 
while experimentally antifcrromagnetioi^, and incorrect 
Fermi surfaces^!. We will thus need to rely on alter- 
nating approaches, such as embedded cluster quantum 
chemical methods, explicitely treating the strong corre- 
lation at the Fermi level. Such methods also present the 
advantages to be able to treat the misfit character of 
the system. For this purpose we will use the same ab 
initio method that gave accurate and reliable results on 
the evaluation of such local effective parameters in other 
strongly correlated systems, that is the embedded frag- 
ment CAS+DDCI methodic. This method yielded very 
good results in systems such as the copper or vanadium 
oxides^, the Na x Co0 2 13 ' 20 , etc., in periodic systems as 
well as misfit ones^i. 

The modulated structures will be taken from refer- 
ence [HI- We will compute the magnetic effective ex- 
change between two nearest neighbors (NN) cobalt ions, 
the effective transfer, the correlation strength and the 
on-site parameters variations that could be responsible 
for electron localization. 



II. METHOD AND COMPUTATIONAL 
DETAILS 

The embedded fragment CAS+DDCI method is a 
spectroscopy method specifically designed to accurately 
treat the strongly correlated character of a set of orbitals 
(designed as active or magnetic orbitals), as well as the 
screening effects on all processes on this set of magnetic 
orbitals. Such a calculation implies large diagonalizations 
and can only be done on a formally finite system. To 
reach the desired goals, the diagonalized configurations 
space includes (I) all configurations that can be built 
within the magnetic orbitals (Complete Active Space or 
CAS), (2) all single and double excitations on each of the 
previous configurations, except for the double excitations 
not contributing to the excitations energies in a pertur- 
bative analysis. Such a choice insures a proper treatment 
of the correlation effect on the cobalt 3d shell as well as 
the screening effects on the cobalt 3d electrons^. 



The embedded fragments are thus composed of a 
quantum part, treated within the CAS+DDCI quantum 
chemical method, and an embedding that reproduces the 
main effects of the rest of the crystal on the quantum 
part. The quantum part should include the cobalt ions 
involved in the sought interaction and their first coordi- 
nation shell. 

The main effects of the rest of the crystal on the quan- 
tum part are double. First, there is the short range ex- 
clusion of the quantum electrons from the space occupied 
by the electrons of the rest of the crystal. This effect 
is modeled using total ions pseudo-potcntialsSS (TIPs) 
that represent the electronic structure of the first shells 
of atoms surrounding the quantum part. The second 
important effect is the Madelung potential. For this pur- 
pose, we used a set of renormalized charges. As far as 
the two periodic directions are concerned, we used the 
real space method developed by A. Gcllc et al££ that 
reproduces the Madelung potential with an exponential 
convergence. For the incommensurate direction, we used 
the same method as in reference l2ll imposing the nullity 
of both charge and dipole moment. 

The question is now "what formal charges should we 
use for the different atoms in order to build the embed- 
ding ?". The calcium and oxygen atoms will be taken 
as Ca 2+ and O 2- ions. The question of the cobalt va- 
lence is a bit more difficult. Indeed, there is no reason 
for the cobalts to present the same valency in the two 
subsystems. Two arguments can help us in finding a 
reasonable evaluation of the charge transfer between the 
two layers. The first argument is a simple comparison be- 
tween the average first neighbor cobalt-oxygen distances 
in the C0O2 layers for CasCo^g and for the Na x Co02 
family. With an average Co-0 distance of 1.904A the 
cobalt valence could be estimated from figure 5 of ref- 
erence Hi] to be similar to the Na2/3Co02 system, that 
is about Co 333+ in the Co0 2 layer, and thus Co 3 05+ 
in the rock-salt layer. The second argument comes from 
a recent evaluation of the cobalt valency using atomic- 
column resolved energy-loss spectroscopy (EELS)2£. In- 
deed Klie et al found a single valency of Co 3+ in the 
rock-salt layers and a mixed valency in the C0O2 lay- 
ers. We thus used a Co 3+ formal charge for the rock-salt 
subsystem and, insuring electro-neutrality, a Co 3 38 for- 
mal charge for the C0O2 subsystem. Such evaluations of 
the cobalts' valences correspond to average values within 
each layer and, at this point, one could wonder about 
the effects of the fluctuations around these values, on 
the Madelung potential acting on the fragments. As a 
matter of example let us evaluate the order of magni- 
tude of the potential modification due to an oxygen va- 
cancy in the rock-salt layer, located in the close vicin- 
ity (~I0A) of the fragment (located in the C0O2 layer). 
The relative modification of the Madelung potential, at 
an atomic position within the quantum fragment, scales 
as SV/V ~ -2/R v /V ~ ICT 5 . The relative fluctuation 
of the potential gradient (the important degree of free- 
dom for our calculations) is even weaker. In addition, the 
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later should be partly compensated by a modification of 
the valency of the cobalt atoms, located in the vicinity 
of the vacancy within the rock-salt layer. One can thus 
expect that such fluctuations, provided their randomness 
and scarcity will little affect the reliability of our results. 

In the quantum calculations, the core electrons where 
treated using effective core potentials^ 6 -, while the valence 
and semi- valence electrons where treated using the asso- 
ciated Gaussian atomic basis set2£. 



related e'^aig 1 and e' g 3 ai g 2 states of the C0O6 embedded 
fragments. Indeed, this excitation energy can be associ- 
ated with the effective orbital splitting of a 3-bands Hub- 
bard model. One sees immediately that whatever the Co 
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III. RESULTS : THE COBALT 3d ORBITALS 

In the C0O2 subsystem of the Ca3Co40g compound, 
the cobalt ions are in a mixed valence state Co 3+ (3d 6 ) 
and Co 4+ (3d 5 ). It is sometimes written in the litera- 
tu re 10 ' 27 that the Co 3+ ions are in a low spin state while 
the Co 4+ ions are in a high spin state. This assumption 
is in contradiction with what was found for the Na x Co02 
family where both cobalt ions are in a low spin state — 
Co 3+ : t 6 g c° and Co 4+ : tLe2- Since these 3d orbitals 
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are responsible for the physical properties, it is of prime 
importance to first understand the cobalt 3d orbital split- 
ting and spin states. 

We thus computed the different spin states of a C0O6 
embedded fragment (Co 3+ as well as Co 4+ ionic states) 
for all crystallographically independent cobalt sites in the 
periodic directions, and a representative set of them in 
the misfit direction. We found that in all cases both the 
Co 3+ and Co 4+ ions are in a low spin state. Indeed, even 
for the Co 4+ valence state the higher spin states are at 
least at 700 mcV (~ 8000 K) above the low spin state. 

As in the Na x Co02 family, the CoOg octahedra are 
compressed along the c axis. One thus expects a split- 
ting of the regular octahedron t 2ff orbitals into a doublet 
e' g and a singlet a\ g . The energetic order between these 
orbitals is crucial since it determines the nature of the 
magnetic orbitals of the Co 4+ ions. A simple crystal 
field evaluation yields the e' g orbitals at the Fermi level. 
However, it was shown that due to the Oh to Sq symme- 
try reduction, the e' g orbitals can and do hybridize with 
the high energy e g ones. It results an energetic stabi- 
lization of the e' g orbitals compared to the a\ g one and 
thus a non degenerated atomic ground state for the Co 4+ 
ions 2 ^ : (e' g ) 4 ai g . In the present system the misfit charac- 
ter of the two subsystems tells us that further symmetry 
breaking occurs and that all cobalt 3d orbitals arc non 
degenerated. 

Figure [2] displays the e aig — e e ' orbital energy split- 
ting as a function of the fourth crystallographic dimen- 
sion parameter (t = raod(nbco0 2 /^RS^)i 1 eN refers 
to the C0O2 subsystem cells in the b direction — see fig- 
ure |3] — &C0O2 an d bus ar e the lattice parameters of the 
two subsystems in the incommensurate direction) associ- 
ated with the incommensurate modulations of the C0O2 
layers. Let us recall that this structural modulation takes 
place along the b direction. These effective orbital split- 
tings where obtained as the energy difference of fully cor- 




Figure 2. (Color online) 3d cobalt orbital energy splittings : 
£a lg — £ e ' g , as a function of the incommensurate modulation 
parameter r. Blue curve with crosses : e aig — £ e ' 1 splitting, 
red curve with plus signs : e aig — £ e ' 2 splitting, black curve 



with stars : average e ai 



splitting. 



site, the a\ g orbital is always of higher energy than the 
e' g ones ; that is the hole is always located in the a\ g or- 
bital. The excitation energy for having a hole located in 
a e'g orbital is always larger than 240 meV ~ 2800 K. The 
average e aig — e e ^ splitting is 261 meV, while the differ- 
ence between the two e' g orbitals ranges between 13 meV 
to 42 meV. Finally one should notice that the variations 
of the e aig — e e i splitting as a function of the different 
crystallographic sites remains small (< ±10meV). 

As mentioned above the relative position of the a\ g 
and e' g orbitals are related to the e' g -e g hybridization. As 
expected from the weak modulation of the orbital energy 
splitting, this hybridization remains stable, close for all 
sites to its average value of 10 degrees. This value is very 
comparable to the values found on different Na x Co02 
systems^. 



IV. RESULTS : ai s -BASED ONE BAND 
MODELS 

The previous section clearly told us that the a\ g -e' g 
orbital splitting is large enough (2800 K) to justify a on- 
band model for our system. In this section, we will thus 
extract from our calculations, both the effective exchange 
integrals, and the parameters of an extended Hubbard 
model. 
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The NN effective exchange can be obtained as the en- 
ergy difference between the singiet and triplet states of 
the C02O10 embedded fragments. The extended Hub- 
bard model parameters necessitate in addition the states 
wave functions information. Indeed, the parameters are 
determined so that the model reproduce both the singlet- 
triplet excitation energy, and the relative weight of the 
dominant configurations, in the states wave functions. 
One should however notice that, not all the parame- 
ters of the extended Hubbard model, can be indepen- 
dently determined from the data. Indeed, only the NN 
hopping integral : t, the difference between the aver- 
age on-site repulsion and the NN repulsion : U — V = 
(Ui + U2) /2 — V12 and the difference between the orbital 
energies and the on-site repulsions : d£ a i g ~ SU can be 
independently determined. 

In this non-periodic system, the models parameters 
were determined for a representative set (sixteen dimcrs 
were computed) of incommensurate distortions values, 
then fitted as a function of the incommensurate dimen- 
sion parameter r. 

r---*-!---_^T- co(i) 



1 / b ' M 
_ , *tt / fC 




I I 
Cellule (n-1) Cellule n 



Figure 3. (Color online) Schematic representation of the dif- 
ferent dimers. 



One can define two types of non equivalent nearest 
neighbor dimers that propagate along the incommensu- 
rate b direction (see figure [3]) ; the A dimers (a — b di- 
rection) and the C dimers (b direction). Indeed, the B 
and D dimers can be deduced from the preceeding ones, 
through the expressions: 



A. Effective exchange and hopping 

Figure Eta) reports the variations of the effective ex- 
change integral J and figure HJb) of the effective hopping 
integral t, as a function of the fourth crystallographic di- 
mension parameter r. Tabic |T] reports the parameters for 
the fit of the computed values according to the formula 

J(t) = J + Ji cos (2tt[t + tp]) + J 2 cos (4tt[t + 
t(r) = t + h cos (2tt[t + £]) + t 2 cos (4tt[t + £]) (3) 



(a) Effective exchange 
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Figure 4. (Color online) (a) Effective exchange, J, and (b) ef- 
fective transfer, t, as a function of the fourth crystallographic 
dimension parameter r. Symbols represent the computed val- 
ues, lines the fits. 



T A ^T B = 0.5-^ - T A (1) 

TC T D = —TC (2) 

In the following discussions, we will refer to "6 direction 
dimers" for the C and D dimers and to "a ± b directions 
dimers" for the A and B ones. 



Let us notice the following points. 

• The NN exchange J may be either antiferromag- 
netic (AFM) or ferromagnetic (FM) (rcsp. negative 
and positive) according to the specific dimer. 

• The NN exchange and hopping arc strongly modu- 
lated for all dimer types. These modulations are es- 
sentially related to the Co-0 distance modulations 
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Dimer 


Jo 


Ji 


J 2 


<P 


to 


h t 2 




A 
C 


2.56 
-1.54 


-6.38 
-0.02 


1.36 
3.49 


0.00 
0.60 


-135.81 
-146.82 


-38.15 12.65 
0.11 27.85 


0.00 
0.10 



Table I. Parameters for effective exchange J and effective 
transfer t fits as a function of the incommensurate dimen- 
sion parameter r. The fitted values are associated with equa- 
tions [3] 



now our computed J values for the CaaC^Og compound 
with both theoretical and experimental results for the 
Na x Co02 family, one sees that the present compound 
- associated with x ~ 0.62 — compares well with the 
Na x Co0 2 family. Indeed, it fluctuating FM/AFM effec- 
tive exchange agrees well with its assumed doping of the 
C0O2 layers close to the phase boundary. 



and not to the Co-Co ones. Indeed, despite the fact 
the Co-Co distances are nearly constant for the C 
and D dimers (b direction dimers), and strongly 
fluctuating for the A and B ones (a ± b directions 
dimers)^, the variation range of both exchanges 
and hoppings present similar amplitudes (^7meV 
and ~13meV respectively). In fact, one should no- 
tice that the amplitude of both the super-exchange 
AFM contribution to J, and the hopping integral, 
decrease when the Co-0 distance increases. 

• J absolute values are very small in amplitude, 
with average values around 2.56 meV for b direction 
dimers, and -1.54 meV for a±6 direction dimers (see 
TableU]). The NN exchange is thus mostly ferromag- 
netic along the b direction and antiferromagnetic 
along the a ± b directions. 



System 


Na x Co0 2 


CaaCoztOg 


X 


0.35 0.5 


~ 0.62 


J (meV) 


-52^ -3622 _n/_i9/_i9/_27^ 


[-2.5,10.3] / [-5.0,2.0] 


Table II. 


Effective exchange values (in 


meV) as a function 



of the C0O2 filing parameter x, present system compared to 
different compounds of the Na x CoC>2 family. 

Comparing these exchange values with the ones ob- 
tained in the Na x Co02 family, one sees that while for 
the superconducting system^ (x = 0.35), the x = 0.5 
and the x = systems^ the exchange are always AFM, 
in the present system both AFM and FM interactions 
were found according to the Co-Co bond. In fact the 
antiferromagnetic character of J decreases with increas- 
ing x (see table |n| for the Na x Co02 compounds. These 
theoretical findings are in agreement with experimen- 
tal data. Indeed, for low x values, the C0O2 lay- 
ers of the Na x CoC>2 systems were found to be anti- 
ferromagnetically coupled^, while for large x values, 
neutrons scattering exhibit A-type antifcrromagnctism^, 
that is ferromagnetic correlations within the C0O2 lay- 
ers which are anti-fcrromagnctically coupled. In fact, 
Lang et al NMR studies^ allowed to propose a magnetic 
phases diagram with a temperature dependant bound- 
ary between the AFM and FM correlations within the 
C0O2 layers. At low temperature the critical doping 
x* is in the 0.6 < x* < 0.7 range and decreases with 
increasing temperature (x*(300K) ~ 0.5). Comparing 



B. Effective on-site repulsion 

FigureOreports the difference between the average on- 
site repulsion on a dimer, U, and the NN repulsion V: 
U - V = (Ui + U 2 ) /2 - Vi 2 . The structural analysis 
shows us that the Co-Co distance variation is very weak 
for the b direction dimers, while it is large for the a ±b 
dimers (See figure [5]). Since the NN repulsion is essen- 
tially related to the Co-Co distances one can assume that 

V is nearly constant along the b direction, but strongly 
varies along the a ± b one. Let us analyze figure O in this 
light. The U — V curve exhibits only weak modulations 
for the C dimers. Since both V and U can be assumed 
as nearly constant for these dimers, the observed modu- 
lations may be related to either of them. In any case one 
can conclude that U remains large and nearly constant 
from site to site (fluctuations range smaller than 0.2 eV 
to be compared to an average U — V of 3.88 eV). The 
sites along the A, B and C, D dimers being the same, it 
results that the strong variation observed on the a ± b 
direction dimers is essentially due to NN repulsion varia- 
tions V(t) associated with the strong fluctuations of the 
Co-Co distances. The amplitude of the variation is quite 
large, with about leV for an average value of 3.5eV (see 
tablc lni] ) . It is not possible using our data to quantify fur- 
ther the NN repulsion, however our results suggest that 

V is large (a few eV in amplitude). An accurate simple 
model should thus take into account the NN repulsion 
between the a± g orbitals and its strong variations ; the 
on-site repulsion being taken as site-independent. The 
average values U — V values found in the present work 
are consistent with the values found for the Na x Co02 
family. Indeed the super-conducting compound exhibits 
U— V = 3.6 eV while it was found for the Nao.sCo02 sys- 
tem 2.6 eV and 2.8 eV according to the crystallographic 
site. 



Dimer 


a 


en 


«2 


c 


A 


3542 


503 


-78 





C 


3878 





-118 


0.1 



Table III. Parameters for U — V fits according to the equation 
U(t) - V(t) = a + aicos(27r[r + C]) + a 2 cos (4tt[t + £]), 
where r is the fourth crystallographic dimension parameter. 
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Figure 5. (Color online) Left scale : difference (meV) between 
the average on-site repulsion on a dimer (U) and the effective 
NN repulsion V as a function of the fourth crystallographic 
dimension parameter r. Right scale : Co-Co NN distances 
(A) as a function of r. 



C. Effective Oi 9 orbital energy 

The on-site repulsion being considered constant, the 
energy difference between NN a\ g orbitals can be ex- 
tracted from our results. Figure [6] reports the effective 
ai s orbital energy : e ai At) (s aig (0) being taken as the 
energy reference). One can see that the amplitude of the 



Effective a lg orbital energy 




Figure 6. ai g orbital energy (meV) as a function of r. The 
r = value is taken as the energy reference. 

e ai (t) modulations is quite large (> 200 meV). Com- 
pared to this modulation, the ai g -e' g energy splitting is 
negligible with an amplitude 10 times weaker. It results 
that the e ai (t) modulation can be seen as an effective 
potential on the whole cobalt atom. Associated with the 
strong NN repulsion V—, one can thus expect a partial 
charge localization, with Co 3+ ions on the low energy 
sites. This result agrees with the observations on the 
Na x Co02 family for which a partial charge localization 
is observed for x > 0.5^1. 



V. DISCUSSION 

Let us now take a closer look on the Seebeck coefficient 
in the Ca3Co40g system. The first point to be stressed is 
the fact that the TEP can be quite different according to 
the intercalation layer (S 300 * (Na 7 Co0 2 ) ^ 100/iV/K 3 ., 
S 300if ([Ca 2 Co0 3 ]o.62[Co0 2 ]) ~ 120 M V/K£io, 
5 300K ([CaBiO 2 ] 2 [CoO 2 ]i. 69 ) ~ 160 M V/K^). The 
second point is the shape of the Seebeck coefficient 
as a function of temperature. At low temperature 
the Seebeck coefficient increases quasi-linearly up to 
~-400K and then saturates to a constant value up to 
about ~400 K 3 - 3 -. These two behaviors can be fitted 
on a correlated Fermi liquid behavior, for the low 
temperature pari— (T < 100 K), and on a Heikes type 
of configurational entropy for the saturating pari— 
(150K< T < 400 K). 

The generalized Heikes formula 3 -^ is a high temperature 
approximation. It supposes that the system is at very 
large temperature compared to the scale of the low energy 
degrees of freedom. In the present compound we showed 
(see section [TTTj) that these degrees of freedom are associ- 
ated with the ai g orbital of the Co0 2 layers cobalt atoms 
3d-shell. Using this formula and the non-degenerated 
character of the magnetic a\ g orbital, one finds (using 
a cobalt valency of 3.38 as discussed previously in sec- 
tion HH, a TEP value of 125 ^V/K. This value is in 
good agreement with the experimental value of 120 /l*V/K 
found for the Ca3Co40g saturation plateau^— One may 
conclude from this agreement that, at room tempera- 
ture, the high temperature limit of the triangular one 
band model based on the a\ g orbitals is already reached. 
However, room temperature only corresponds to 25 meV, 
while the band width of a tight binding model on a tri- 
angular lattice is W = 9|i| ; using our present transfer 
estimations (\t\ > 100 meV) it comes a band width of the 
order of 10000 K. How can we account for such a large 
difference in order of magnitude? Why does the Heikes 
model fits so nicely the room temperature Seebeck coef- 
ficient for the Co0 2 -based compounds? 

Different arguments can be advanced in order to ex- 
plain this high temperature behavior of the TEP at 
300 K. 

Let us first analyze the behavior of the Seebeck coeffi- 
cient, as a function of temperature within a one-band 
tight binding model on a triangular lattice. The Boltz- 
mann equation yields the following Seebeck expression 



S(T) 



k B T 



de at 



Je^de 



k B T J ^de 
S tr {T) + S^T) 



(4) 

(5) 

(G) 
(7) 



where fo(T, e) is the Fermi function, e the orbitals energy, 
fi the electrochemical potential ; the collision time as well 
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as the degeneracy is supposed energy independent. Fig- 
ure [7] reports S(T) and its two components St r {T) and 
Sfj,(T). At low temperature St r (T) dominates, reaches a 
maximum for T/W ~ 1/4, then slowly decreases toward 
zero in the high temperature limit. On the contrary the 
cntropic part S^fT) increases for all temperatures and 
reaches its high temperature limit for temperatures of 
the order of a few band widths. The sum of the two 
terms results in a Scebcck coefficient reaching its satura- 
tion value at much lower temperature than its entropic 
part. Indeed for T/W = 0.5, S is already larger than 
90% of its high temperature value. This simple analysis 
already allows us to reduce, from a few band widths to 
half the band width, the estimation of the temperature 
at which the TEP saturation should be observed. This 
is however not enough to account for the observed satu- 
ration temperature. Indeed W/2 = 9\t\/2 still accounts 
for about 5000 K. 

A second argument is the extremely large value of the 




T/W 

Figure 7. (Color online) Seebeck coefficient of a triangular 
lattice as a function of the temperature, normalized to the 
band width. Blue dashed curve : the entropic term : S^(T), 
red dot-dashed curve : the Fermi liquid part : Str(T). 

correlation strength in these compounds. Indeed, our cal- 
culations yield for the average U/t a value of 26. Despite 
this very large U/t, transport measurement ;) 6 ' 10 find the 
system metallic. One thus expect its Fermi level quasi- 
particle pic to be strongly renormalized due to the strong 
correlation effects^. 

Finally let us remember that our results suggest a par- 
tial electron localization due the large site energy mod- 
ulation (> 200 meV) and a large value of NN repulsion. 
It results that the lattice pattern seen by the Fermi level 
electrons can be expected to be quite different from the 
homogeneous triangular lattice assumed in these systems. 
This effect is also expected to strongly reduce the quasi- 
particlc bandwidth at the Fermi level. 
All together, it seems quite possible that the band- 
width reduction is large enough (20 to 30 times) so that 
0.5W o s ~ 200 K, and thus that the high temperature 
Heikes model appears valid for the present (and related) 
systems. 



Let us now take a look on the TEP behavior at tem- 
peratures larger than 400 K. Indeed, at T c = 396 K the 
Ca3Co40g compound experiences a first-order metal- 
semiconductor phase transition associated with a linear 
increase of the Seebeck coefficient for T > Tj&. Dif- 
ferent hypotheses have been proposed in order to ac- 
count for this phase transition and the high temperature 
phase properties. The first one associates the transition 
with a cobalt spin-state transition^ -. However our cal- 
culations show that the first atomic excited spin state 
of the cobalt atom is at least at 700 meV, i.e. 8100 K, 
above the ground state. Such an energy difference is 
clearly incompatible with a phase transition energy scale 
of 400 K. The second hypothesis associates the transi- 
tion with a change in the cobalt t2 g orbitals occupancies. 
However this hypothesis does not resist to the energy 
scale of the a\ g -e' g energy splitting found in section IIIII 
Indeed, the 240meV~2800K is still far too large for the 
400 K energy scale of the transition. The last hypothe- 
sis is an electron localization of the e' g orbitals, associ- 
ated with a structural phase transition^. Even if the 
assumption that the e' g orbitals are the Fermi level ones 
has been proven incorrect (see reference [H or for the 
present system section iHTj) . the idea of an electron local- 
ization induced by an increase of the in-plane distances 
in the C0O2 layers seem quite possible and would explain 
the metal-semi-conductor phase transition. Muguera and 
Grebille^ solved the crystallographic structure in the 
high temperature phase. They showed that the tran- 
sition is due to a rearrangement of the central CoO layer 
of the Ca2Co03 subsystem. This rearrangement acts 
on the C0O2 layers through the O-Ca interaction be- 
tween the two subsystems. It results in a modification 
of both the cobalts and oxygens modulations (but not 
of the average structure) along the a and c directions. 
While such a weak structural modification can account 
for the metal-semiconductor transition (due to the large 
correlation strength and associated closeness to the Mott 
transition), it cannot account for the associated increase 
of the Seebeck coefficient. Indeed, such an increase ne- 
cessitate the opening of new degrees of freedom at the 
Fermi level and the weakness of the structural modifica- 
tions within the C0O2 layers cannot account the opening 
of either the higher cobalt spin states or the e' g orbitals. 
This study however show us that there is another degree 
of freedom up to now ignored that should be taken into 
account, that is the partial charge localization due to the 
action of the subsystem modulation over the Fermi level 
orbital energies. Indeed a modification of the modula- 
tion vector will be expected to modify the energy profile 
pictured in figure [6l thus the charge localization configu- 
ration and the lattice pattern seen by the Fermi electrons 
at the given temperature. Out of any other energetically 
compatible hypothesis, the latter seems the most plausi- 
ble. 

Let us now extend our discussion to the other impor- 
tant parameters for the thermoelectric performance. As 
we all know, a good figure of merit is associated with a 
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large Seebeck coefficient S, a low resistivity p and a low 
thermal conductivity k. 

S : We saw in the preceeding paragraphs that the 
Hcikes formula is fully justified to describe the See- 
beck coefficient of the layered cobalt oxides (pro- 
vided a large U/t value and a partial charge lo- 
calization). The Heikes formula is a pure entropic 
one and only depends on the band filling of the a\ g 
cobalts magnetic orbitals of the C0O2 layer. Thus, 
as far as S is concerned, the misfit character es- 
sentially acts (i) on the filling of the C0O2 layer, 
since the latter should adjust to the incommen- 
surate ratio (bco0 2 /bRs) in order to preserve the 
electro-neutrality of the system, (ii) on the renor- 
malization of the a\ g band width (through a par- 
tial localization). In both cases these are indirect 
effects (acting through the electrostatic interaction 
between the two layers) and the effect of the spe- 
cific modulation vectors, or even their incommen- 
surate character, will be expected to modulate the 
S value, but not its order of magnitude. 

p : In these Co02-based systems, the observed low re- 
sistivity is essentially due to the triangular pat- 
tern of the transfer interactions in the C0O2 layer. 
Indeed, in a triangular lattice, the Mott tran- 
sition is repelled at very high correlation val- 
ues, and the system remains essentially metal- 
lic (as experimentally observed in the misfit 
cobalt oxides which present low metal to insu- 
lator transition temperatures, between 50K-80K 
for the present compound^ ' 10 ' 40 , around 20K for 
the [Sr 2 Coi_ x Bi x 03][Co02]i.8 compounds^ or the 
[Bii. 95 Bai.95Rho.i04][Rh02]i.s^, or even remain- 
ing metallic at all temperatures such as the 
[Sr 2 02]o.53[Co0 2 )^). As a matter of illustration, 
let us compare the U/t ~ 1522, ratio found for the 
metal to insulator transition of a half-filled, trian- 
gular, Hubbard model, and the U/t = ratio for 
the same transition when the transfer interactions 
pattern is square. The misfit character can thus be 
expected to be non relevant compared to the basic 
triangular interaction pattern as far as the resis- 
tivity is concerned. At a higher order of precision, 
however, it is expected to sightly increases the re- 
sistivity, through the induced partial localization. 

k : The main effect of the misfit character is most 
probably on the thermal conductivity. Indeed, 
one can expect that the non periodicity and the 
associated modulated potential seen by otherwise 
equivalent ions, induce a partial localization of 
the phonons dispersion and thus a lowering of 
the thermal conductivity. This point is however 
difficult to quantify since it would necessitate a 
comparative study of the phonons dispersion spec- 



tra for different misfit and non misfit Co02-bascd 
compounds. 

VI. CONCLUSION 

We performed ab initio calculations of the effective 
local parameters of the thermoelectric Ca3Co40g com- 
pound, taking into account the incommensurate distor- 
tions. As in the Na x CoC>2 family, cobalts ai g -e' g orbitals 
splitting is large so that the low energy physics can be 
accounted for using the a\ g orbitals only. We showed 
that the incommensurate modulation of the C0O2 layer 
cannot be neglected in order to understand the com- 
pound physical properties. Indeed, site energies, NN ex- 
change, NN hopping as well as first neighbor repulsions 
are strongly modulated due to the misfit character of the 
system. It particular the site energies large modulation 
let us expect a partial Fermi electrons localization. The 
effective exchange remains weak in amplitude, however 
fluctuates between antifcrromagnctic coupling and ferro- 
magnetic one. This results can be compared with the 
on-plane AFM-FM transition line in the Na x Co02 fam- 
ily associated at T = K with similar cobalt valency as 
found in the present compound. 

We discussed the temperature behavior of the See- 
beck coefficient and showed that the cobalts a\ g orbitals 
bandwidth is strongly renormalized, due to very strong 
electron correlation (U/t ~ 26), and the modification of 
the effective lattice pattern seen by the itinerant elec- 
trons, associated with the partial charge localization. 
This strong bandwidth renormalization is most proba- 
bly responsible for the plateau due to the high tempera- 
ture behavior observed in the Seebeck coefficient on the 
[150 K, 400 K] range. We further showed that he sudden 
increase of the TEP for T > T c = 396 K, cannot be at- 
tributed either to a cobalt spin state change or to a mod- 
ification of the dig vs e' g orbitals occupation. In fact the 
396 K metal-semiconductor transition is associated with 
a structural rearrangement within the Ca2CoC>3 subsys- 
tem inducing a modulation modification within the C0O2 
layer. This modification induces a site energy pattern 
change, and thus a localization pattern change, which 
we think to be the best candidate for the change in the 
Seebeck behavior. 
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